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ABSTRACT 

We show that the prompt and afterglow X-ray emission of GRB060218, as well as its early (t < 1 d) 
optical-UV emission, can be explained by a model in which a radiation- mediated shock propagates 
through a compact progenitor star into a dense wind. The prompt thermal X-ray emission is produced 
in this model as the mildly relativistic shock, (3 w 0.85 carrying few x 10 49 erg, reaches the wind 
(Thomson) photosphere, where the post-shock thermal radiation is released and the shock becomes 
collisionless. Adopting this interpretation of the thermal X-ray emission, a subsequent X-ray afterglow 
is predicted, due to synchrotron emission and inverse-Compton scattering of SN UV photons by 
electrons accelerated in the collisionless shock. Early optical-UV emission is also predicted, due to the 
cooling of the outer SM ~ 10 _3 M Q envelope of the star, which was heated to high temperature during 
shock passage. The observed X-ray afterglow and the early optical-UV emission are both consistent 
with those expected in this model. Detailed analysis of the early optical-UV emission may provide 
detailed constraints on the density distribution near the stellar surface. 
Subject headings: gamma rays: bursts — supernovae: general -shock waves 



1. INTRODUCTION 

In our previous paper (jCampana et al. 1 [2006) the 
points mentioned in the abstract were outlined only 
briefly, using order-of-magnitude arguments and with 
very little explanation, due to space limitations. Here 
we present a more detailed explanation and analysis of 
the prompt thermal X-ray emission and X-ray afterglow, 
as well as a calculation of the early optical-UV emission, 
and show that some claims made in recent publica t ions 
dChisellini et al. Il200l iFan et al. I l2fM Eil l200fih . ac- 
cording to which the observations are inconsistent with 
the massive wind interpretation, are not valid. 

GRB060218 was unique mainly in two respects: it 
showed a strong thermal X-ray emission accompanying 
the prompt non-thermal emission, and a strong optical- 
UV emission at early, t < 1 d, time. We show here 
that these features, as well as the X-ray afterglow, can 
all be explained by a model in which a radiation medi- 
ated shock propagates through a compact progenitor star 
into a massive wind . We have shown in a separate paper 
(|Wang et al. II2006T ) that the prompt non-thermal X-ray 
emission can also be explained by this model. 

As detailed below, the prompt thermal X-ray emis- 
sion can be explained as shock breakout at a radius of 
~5x 10 12 cm, which requires a Thomson optical depth 
(of the plasma ahead of the shock) of r w 1. Breakout 
may occur at ~ 5 x 10 12 cm if this is the stellar pro- 
genitor radius. However, since the progenitor is likely to 
be smaller, we suggested the possibility of it being sur- 
rounded by an optically thick wind. Another possibility 
is a pre-explosion ejection of a small mass, ~ 10 _6 Mq, 
shell. Here we shall adopt the wind interpretation, since 
it allows one to explain also the X-ray afterglow. So far, 
no other quantitative physical models have been worked 
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out for t he thermal X-ray emissi on, n or for the X-ray af - 
terglow. ISoderberg et all (|2006l ) and IFan et al. I (|2006f ) 
have suggested that the afterglow X-ray emission is due 
to extended activity of the source, for which th ere is no 
model or explanation. iGhisellini et al. I (|2006l ) suggest, 
for explaining the X-ray afterglow, an ansatz consisting 
of the ad-hoc existence of electrons with some prescribed 
energy distribution, which is different at different times 
to account for the observations, without a model for the 
dynamics of the plasma or for the electron energy distri- 
bution. 

We not e that the radio a ftergl ow of GRB060218 dis- 
cussed by ISoderberg et aT] (|2006| ) is difficult to explain 
within the co ntext of the current m odel. Indeed, as 
pointe d out bv lSoderberg et all (2006) and bv lFan et al. I 
(2006), it is difficult to explain the radio afterglow and 
the X-ray afterglow of GRB060218 as due to emission 
from a single shock wave. These authors have thus 
chosen to construct models that account for the radio 
emission only, attributing the X-ray afterglow to a con- 
tinued activity of the source of an unexplained nature, 
and not accounting for the prompt X-ray emission and 
for the early optical-UV emission. We adopt a differ- 
ent approach, showing that the prompt (thermal and 
non-thermal) X-ray emission, the early optical-UV emis- 
sion, and the late X-ray afterglow can all be explained 
within the context of the same model. We argue that 
it is the radio afterglow, rather than all the other com- 
ponents, which remains unexplained, and which should 
be attributed to a different component. Since the radio 
emission carries only a negligible fraction of the energy, 
and given the large anisotropy of the explosion, it is not 
difficult to imagine the existence of such an additional 
low energy component. 

In § [2] we discuss the observations and our model. In 
§[3] we compare our analysis with those of other authors. 
Our results and conclusions are discussed in § |4j 

2. OBSERVATIONS AND MODEL 

2.1. Thermal X-ray Emission: Shock Breakout 
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Possibly the most distinguishing feature of GRB060218 
is the strong thermal X-ray emission accompanying the 
prompt non-thermal emission. The temperature of the 
thermal X -ray photons observed up to 3 x 10 3 s is T w 
0.17 keV (|Campana et al. 112006ft . The integrated flux 
of the black-body X-ray component, 6 x 10 _6 erg/cm 2 
(see figured]), corresponds (using d = 145 Mpc and after 
correcting for the flux that falls outside the XRT band) 
to a thermal X-ray energy of E t h ~ 2 x 10 49 erg. This 
is only an approximate estimate, and the total thermal 
energy may be somewhat larger, due to the gap in XRT 
observations between 3 x 10 3 s and 6 x 10 3 s. We consider 
here a model where this radiation is due to a "shock 
breakout" . 
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Fig. 1. — Upper panel: the X-ray light curve (0.3-10 keV). Open 
black circles mark Swift-XRT data, rilled squares Chandra data 
(Soderberg et al. 2006) and filled circle XMM-Newton data (de 
Luca & Marshall 2006). Open squares denote the black body com- 
ponent. A smooth-burst with exponential decay fits the prompt 
part (solid line). A power law decay with index -1.25 well de- 
scribes the afterglow tail (dashed line). Count rate-to-flux conver- 
sion factors were derived from time-dependent spectral analysis. 
The X-ray light curve has a long, slow power-law rise followed 
by an exponential (or steep power-law) decay. At about 10,000 s 
the light curve breaks to a shallower power-law decay with index 
— 1.2 ± 0.1. Lower panel: the UVOT light curve. Filled circles of 
different colors (on-line version) represent different UVOT filters 
(from bottom to top between 104 — 105 s): V (centered at 544 
nm); B (439 nm), U (345 nm), UVW1 (251 nm); UVM1 (217 nm) 
and UVW2 (188 nm). Specific fluxes have been multiplied by their 
FWHM widths (75, 98, 88, 70, 51 and 76 nm, respectively). 



Supernova (SN) shock waves become radiation- 
med iated when propagating th rough the stellar envelope 
(see IWoosley fe Weaver Ill986l for review). The energy 
density behind the shock is dominated by radiation, and 
the mechanism that converts kinetic energy to thermal 
energy at the shock transition is Compton scattering. 
The optical depth of the shock transition layer is r s ~ 
c/v s , where v s is the shock velocity. This thickness is de- 
termined by the requirement that the time it takes a fluid 
element to flow through the transition layer, ~ A s /v s , 
should be comparable to the diffusion time of photons 
across this layer, ~ t s A s /c. As the shock propagates 
outward, the (Thomson) optical depth of the plasma ly- 
ing ahead of the shock decreases, and when this optical 
depth becomes comparable to t s , Compton scattering 
can no longer maintain the shock. At this point, radi- 



ation escapes a head of the shock, producing the shock 
breakout flash ([Colgate 111974 iKlein fc Chevalier1ll978t 
lEnsman fe Burrows I Il992t iMatzner fe McKee I 1 1999ft . 
and the shock becomes collisional (or collisionless, see 
IWaxman fe Loeb Il200lh . 

Hereafter we use the term "shock breakout" to denote 
the event of the transition from radiation to collisional 
(collisionless) shock mediation, accompanied by the emis- 
sion of radiation. If the optical depth of the wind sur- 
rounding the progenitor star is small, shock breakout will 
take place as the shock approaches the stellar surface. It 
is for this reason customary to identify shock breakout 
with the emergence of the shock from the stellar surface. 
However, if the optical depth of the wind is large, > r s , 
breakout would occur once the shock reaches a radius 
where the wind optical depth drops to r s . 

As argued in ICampana et alH (|2006l ). in order to ob- 
tain a breakout flash with energy of E t h — 10 49 erg and 
temperature T ss 0.17 keV, breakout must occur at a 
radius of r ~ 5 x 10 12 cm, and the shock must be mildly 
relativistic (which implies r s ~ 1). Since the progeni- 
tor star is presumably smaller than this (e.g. if it is a 
Wolf-Rayet star), a value of r ~ 1 at r ~ 5 x 10 12 cm 
may be obtained either by assuming that the progenitor 
is surrounded by a dense wind, or that an outer shell of 
the star was ejected to this radius prior to the GRB ex- 
plosion. The mass of the shell required to obtain t s ~ 1 
is only ~ Aitr 2 / k ~ 1O -6 M0, where K = 0.2g/cm 2 is 
the Thomson opacity for ionized He. For the calculation 
below we adopt a density profile p oc r~ 2 , as would be 
expected for the wind model. The results obtained for 
the breakout radius, shock velocity and plasma density 
at this radius are not sensitive, however, to the details of 
the density profile shape. 

We consider therefore a shock wave of velocity (3c 
driven into the wind surrounding the progenitor, whose 
4- velocity is u = 7/? where 7 = (1— (3 2 )~ x l 2 is the Lorentz 
factor. Let us first derive the energy and temperature of 
the post-shock radiation. The post-shock temperature 
Td is related to the post-shock pressure by aTj = 3p, 
and the observed temperature is given by "fdTd, where 
jd = (1 — fid)~~ 1 ^ 2 an d ftdC is the post-shock plasma ve- 
locity {(3d < (3). For a strong radiation dominated shock 
p = fu 2 pc 2 , where p is the pre-shock density and f(u) 
is a factor of order unity. As we show below, the shock 
is required to be mildly relativistic, u — a few, for which 
we may approximate / = 0.8 (for h>1, / = 2/3), and 
(u d /u) 2 w 0.6 ((u d /u) 2 = 1/2 for u > 1). This defines 
the observed temperature, T, to be 



T~l U 3 / 2 (a-V 2 ) 1/4 



(1) 



(for u 1, the numerical coefficient is 1/2 1 / 4 instead 
of 1). The energy carried by the radiation may be esti- 
mated by noting that the energy density of the radiation 
is given, in the observer frame, by — p + 7 2 (aT' ( | + p) = 
(4u^ + 3)p, and that the thickness of the shocked plasma 
shell may be approximated as d = {(3 — (3d)r, where r is 
the shock radius. This yields 
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(for u 3> 1, the numerical coefficient is 2/3 instead of 
0.5). 
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For a wind density profile, the optical depth of the 
wind at radius r is simply r(r) = np(r)r. Expressing the 
density in terms of r(r) and using eqs. |T]) and (2|) we 
find that the shock velocity at breakout is 
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and that breakout occurs at a radius 
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where £th = 10 i?th,49 erg and we have used k = 
0.2g/cm 2 and r = 1 in the numerical evaluations. Note 
that the results depends only weakly on the exact values 
of k and t. 

Eqs. (3]) and (HJ) imply that the breakout of the shock 
takes place at a radius ~ 7.8 x 10 12 cm, and that the 
shock is mildly relativistic at breakout, u 2 — 2.5 i.e. 
/3 = 0.85. These results are in agreement with the orde r 
of magnitude estimates given in lCampana et al~l (2006). 
Since the shock is found to be mildly relativistic, using 
r = 1 at shock breakout is justified. In order to produce 
the shock that is driven into the wind, the explosion is 
thus required to produce a mildly relativistic shell, [3 = 
0.8, with an energy of a few x 10 49 erg. This is remark- 
ably similar to the case of GRB980425/SN1998bw, for 
which the ejection of a shell of energy 10 49 7 erg and ve- 
locity (3 == 0.8 was inferred from X-ray dWaxman Il2004bf) 
and radio feulkarni et al. lll998t | Waxman fc Loeb Ill99£ ; 
IChevalier fc Li 1119991 : IWaxman 2004a) observations. 

The wind density at the breakout radius, p === 1/nr = 
10~ 12 g/cm 3 , is a lso in agreement with our results in 
ICampana et al. I (2006). It is straightforward to ver- 
ify that the energy density is dominated by radiation, 
aT A /nT ~ 10 9 . For a wind velocity of ~ 10 3 km/s it cor- 
responds to a mass loss rate of few xl0 _4 MQ/yr. The 
relevant mass loss here is that occurring within a day 
or less of the explosion. The data currently available on 
wind mass losses suggesting M < few x 10 _4 MQ/yr refer 
to m uch longer timescales, for s tars well before any explo- 
sion (|Mevnet fc Maeder Il2007f ). Physically, it is however 
quite plausible that the mass loss increases considerably 
as the evolution of the core rapidly approaches the final 
collapse, accompanied by a rapid increase in the lumi- 
nosity and the envelope expansion rate. 

Note that if both the shell and the wind were spher- 
ically symmetric, the characteristic timescale would be 
Rph/c = 260 s (for the inferred post shock velocity, 
IdPd = 1-2, the effects of relativistic beaming are not 
significant), while the obs erved timescale of the t hermal 
X-ray emission is ~ 10 3 s (jCampana et al. II2006D . How- 
ever, an anisotropic shell ejection is a natural expecta- 
tion in a core collapse GRB, si nce strong rotation is a 
requis ite to make the jet (e.g. [MacFad ven fc Wooslev I 
(1999)). Even "normal" core collapse SN simulations 
show strong rotation-related anisotr opy in the expand- 
ing gas (e.g. iBurrows. et al. I (|2007l )). Thus, the semi- 
relativistic outer shell ejected is likely to be anisotropic, 



either due to an anisotropic explosion, or due to being 
driven by a jet. This, a s well as an anisotropic w ind pro- 
file caused by rotation (Mc vnet fc Maeder~ll2007l ) should 
lead to significant departures from sphericity in the shock 
propagation. Anisotropy is in fact a prediction of this 
model, whi ch is supported by the d etection of linear po- 
larization (|Gorosabel et al. 1I2006D . In an anisotropic 
shock, however, the timescale is no longer the naive 
spherical r/c value, but is rather given by the sideways 
pattern expansion timescale, which depends on the angu- 
lar velocity profile of the anisotropic shell (e.g. at larger 
angles the shock emerges later due a decreasing velocity 
profile or due to an increasing wind density away from 
the symmetry axis, etc.). 

2.2. X-ray afterglow: Wind-shell interaction 

If the wind shock breakout interpretation is adopted, 
the subsequent interaction of the ejected relativistic shell 
with the wind is expec ted to produce an X-ray after- 
glow (W axmanl [2004b) . The electrons accelerated to 
high energy in the collisionless shock driven into the wind 
emit X-ray synchrotron radiation, and may also inverse- 
Compton scatter optical-UV SN photons to the X-ray 
band. 

As the shock wave driven into the wind expands, it 
heats an increasing amount of mass to a high tempera- 
ture and the initial energy of the ejected shell is trans- 
ferred to the shocked wind. The shell begins to decelerate 
beyond a radius at which the shocked wind energy be- 
comes co mparable to the ini tial shell energy. This occurs 
at a time (jWaxman Il2004bft 



El 49 



(5) 



Here, Ek == 10 49 EkA9 erg is the kinetic energy of the 
shell and 



M/v v 



10- 5 (M Q /yr)/10 3 (km/s) 



(6) 



For the wind density and shell velocity inferred in § 12.11 
P = 0.85 and m > 10, td ec ^ 10 4 s. At t > td ec the energy 
is carried by shocked wind plasma, and is continuously 
transferred to a larger mass of newly shocked parts of the 
wind. Thus, the rate at which energy is transferred to ac- 
celerated electrons is ~ e e Ek/t, where e e is the fraction 
of the post-shock thermal energy carried by electrons. 
If the electrons cool on a time scale shorter than the 
expansion time, ~ t, this would lead to a bolometric 
luminosity Lb ~ 10 43 e ej _ii?fc j 49(t/lday) _1 erg/cm 2 /s. 
The fraction of this luminosity emitted as X-rays de- 
pends on the electron energy distribution. For elec- 
trons accelerated to a power-law distribution in energy, 
c?logn e /dlog7 e w -2, vh v w L B /2 log{~f max / j min ) w 
l0 A2 e e - 1 E kA9 (t/lday)~ 1 erg/cm 2 /s, where ^ max and 
Jmin ar e the maximum and minimum electr on Lorentz 
factor. The more detailed analysis given in Waxman 
(2004b) of the s ynchrotron emiss ion of shock accelerated 
electrons gives (jWaxman Il2004bl eqs. (6) and (7)) 



for 



fijAsynch ~ 3 x 10 e e -iE kA9 t d erg/s (7) 
v > v c re 3 x 10 10 e B 3 Z 2 1 (77i/ 1 0) _3/2 ^Hz. (8) 
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Here, t — ltd day, e e = 10 e e ,-i and €b = 10 es,-i 
are the fractions of the post-shock thermal energy carried 
by electrons and magnetic field and v c is the cooling fre- 
quency, the frequency of synchrotron photons emitted by 
electrons for which the synchrotron cooling time equals t 
(higher energy electrons cool faster and emit higher en- 
ergy photons). Eq. ([7]) implies an X-ray luminosity in 
Swift's XRT band, 0.3-10 keV, of 



J X, synch 



io 42 £r 



L-Efc,49*d ler S/ s - 



(9) 



Let us consider next the inverse- Compton scattering 
of SN optical-UV photons by the shock-heated electrons. 
The inverse- Compton luminosity emitted by electrons 
with Lorentz factor 7 e is given by Lie ~ lt T {le)LsN , 
where r(7 e ) is the Thomson optical depth of these elec- 
trons and Lsn is the SN luminosity. The lowest energy 
IC photons are produced by the lowest energy, ther- 
mal, electrons, {hv)ic,T ~ (47^/3)2 eV, where jt is 
the Lorentz factor of the lowest energy electrons. 7t 
is determined from the post-shock thermal energy den- 
sity, [2/(r 2 — l)]pv 2 s» pv 2 where v s is the shock velocity 
and p — M /4nr 2 v w (F w 5/3 is the adiabatic index), 
through An e 7TTO e c 2 = e e pv 2 - Here, A is a dimensionless 
constant, the value of which depends on the exact form 
of the electron energy distribution. For a power-law dis- 
tribution, dlogn e /dlog7 e w -2, A « log(7 maa; /7 mi „). 
The shock velocity ma y be inferred from the shock ra- 
dius, which is given by (jWaxman Il2004bl eq. 2) 
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Using v s = (2/3)R s /t we find 
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Thus, on a time scale of a days IC scattering by the low- 
est energy (thermal) electrons is expected to contribute 
to the X-ray flux. It is important to note here that the 
electrons producing the X-ray synchrotron flux, eq. 
are of much higher energy, 7 e 3> 7t, and lie at the high 
energy part of the accelerated electron energy distribu- 
tion. 

The Thomson optical depth of the thermal electrons is 
approximately given by kM (r) / Airr 2 = npr, where M(r) 
is the wind mass accumulated up to r. Assuming that the 
thermal electrons do not lose all their energy by IC scat- 
tering on a time scale shorter than the expansion time t, 
the IC luminosity produced by the thermal electrons is 



Lic,t ~ Kprj T L S N 
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The IC luminosity given by eq. (|13p may exceed, depend- 
ing on the values of A and e e , the luminosity given by 



eq. ©. Ljc.t > Lx, synch is not, of course, a valid result, 
since L v, S ynch given in eq. |9|) is the luminosity obtained 
assuming that the electrons lose to radiation all the en- 
ergy they gained from the shock. The IC luminosity is 
thus limited by Lx.synch, and L ic ,t > L x , sy nch simply 
implies that the thermal electrons lose all their energy to 
IC scattering on a short time scale. 

Let us compare the predicted X-ray afterglow with 
the observed one. The observed X-ray afterglow, 
following the prompt emission which ends at ~ 
10 4 s, is well approximated by (see figure [I]) fx = 
10 _12 t ti " 1 ' 2±0 ' 1 erg/cm 2 /s, which corresponds to Lx = 
2 x 10 42 <7 1 ' 2±01 crg/s. This is in excellent agree- 
ment with the predictions of eqs. (J9j) and (fT3|) for 
Ek = few xlO 49 erg. The fact that the energy of the 
shock driven into the wind inferred from the X-ray after- 
glow, Ek, is comparable to that of the thermal X-rays, 
E t h, supports our model, in which both are due to the 
same shock driven into the wind. 

At early time, t < 1 day, the emission is ex- 
pected to be dominated by the synchrotron component, 
and for the power-law energy distribution of electrons, 
o?logn e /dlog7 e w —2, the X-ray spectrum is expected 
to follow vL v cx u°. On a time scale of a few days, the 
emission is expected to be dominated by IC scattering of 
thermal SN photons by thermal shock electrons. At this 
stage, the X-ray spectrum may become steeper, reflect- 
ing the energy distribution of the lowest energy electrons 
heated by the shock. These results are consistent with 
the vL v cx v X-ray spectrum measured at t ~ 1 day, 
and with the indication, based on an XMM-Newton ob- 
servation, that at a later ti me, t ~ 3 day, the spectrum 
is steeper, vL v cx z,-i-3±o.6 (jde Luca fc Marshall Il2006| ). 

The following point should be made here. As the shock 
speed R s decreases with time, shells ejected with veloci- 
ties lower than that of the fast, (3 ~ 0.8, shell may "catch 

— 1/3 

up" with the shock [R s /ct — 0.6i d for £^,49 = 1 and 
m = 10). This may lead to increase with time of the 
shock kinetic energy Ek, and thus to a modification of 
the X-ray light curve. The fact that the X-ray luminos- 
ity decreases roughly as 1/i up to ~ 10 d implies that 
the energy in shells ejected with (3 > 0.3 is not much 
larger than that of the fast shell, i.e. not much larger 
than 10 49 erg. 

2.3. Early optical-UV emission: Envelope cooling 

As the SN shock propagates through the stellar enve- 
lope, it heats it to ~ 1 keV. As the envelope expands, 
the photosphere propagates into the envelope, and we 
see deeper shells with lower temperatures. We derive 
here a simple analytic model for the photospheric radius 
and temperature, based on which we can derive approx- 
imately the flux and temperature of the escaping radia- 
tion. We assume that the density of the stellar envelope 
near the stellar surface is given by 



Po(r) = P1/2S" 



(14) 



where S = (1 — r/R*) (i?* is the stellar radius), and 
that the the plasma ahead of the photosphere is highly 
ionized He, so that k — 0.17cm 2 /g. The velocity 
of the shock as it propaga tes through the envelope 
is approximately given by (jMatzner fc McKee I [1999; 
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iTan. Matzner fc McKee1l2001h 



E, 



■0.8 



E, 



1/2 



1/2 



M, 



0.2 



pr° 



5 



-0.2ri 



(15) 



where M e j is the mass of the ejected envelope, and E e j 
is the energy deposited in the envelope. Since the shock 
is radiation dominated, the energy density behind the 
shock is given by u rad ,o/3 = {6/7)p v 2 , 

1 8 

aT 4 {5) = u iad! o(§) = yPo(S)v s (5) 2 . (16) 

After the shock breaks through the envelope, the enve- 
lope begins to expand and cool. It is useful do label the 
shells with Lagrangian coordinates, defining S m (S)M e j 
as the (time independent) mass that lies ahead of a shell 
origin ally located at r = (1 — 8)R*. Mat zner fc McKeel 
( 1999) show that the final velocity (after acceleration due 
to the adiabatic expansion) of each shocked shell, Vf, is 
related to the shock velocity approximately by Vf(5 m ) = 
2v s [S(S m )]. After significant expansion, Vft ^> the 
radius of each shell is given by r(S rn , t) — Vf(5 m )t. Given 
r(8 m , t) it is straightforward to derive the time dependent 
shell density. The shell's (time dependent) temperature 
is then determined by T/Tq — (p/po) 4 ^ 3 , which holds for 
adiabatic expansion. After some tedious algebra we find 
that, for Vft ^> the photosphere is located at 



<W(*) = 3.8 x 10" 3 /- 07 



E o.s 



(M ej /M Q ) 



(17) 



Here, E e j = 10 5 E^si erg and f p = P1/2/P0 is the ratio 
between pi/ 2 and the average envelope density p~Q. Al- 
though f p depends on the structure of the progenitor star 
far from the surface, where eq. (fT4|) no longer holds, the 
results are very insensitive to its value. Using eqs. (115[) 
we find that the radius of the photosphere is 



r ph (t)=3.2xl0 1 V° 4 



F 0A 



(M ej /M Q )™ 



t , 8 cm, 



(18) 



and using eq. (|16[) we find that the temperature of the 
photosphere is 



T ph (t) = 2.2/; 



0.02 



77.0.02 



(Mej /Mq ) 

H ere, R* = 10 12 -B* ,i 2 cm 
ICampana et al. I (|2006l) find R 



pl/4 ,-0.5 „v 
0.03 U *A2 T d eV - 



(19) 



3 x 10 14 cm and 



T « 3 eV at t = 10 s. This is clearly consistent with 
the cooling envelope interpretation. We note that at 
t = 10 5 s the emission originates from a shell of mass ~ 
10 ~ 3 MfT), in excellent agre ement with our rough estimate 
in (jCampana et al. 1l2006l ). At the time t <C 10 5 s opti- 
cal photons are in the Rayleigh- Jeans tail, and the flux 
is approximately given by f v ~ 2ir(i//c) 2 T p Y 1 (r p \ 1 / D) 2 , 
where D is the distance to GRB060218. Using eqs. (fTgj) 
and (fT9l) we find 
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vf v « 1.3 x 10- iJ / p 



-^,51 R l/4 
0.6-"*, 12 



(M ej /M Q ) 



5 x 10 14 Hz 



t\ A erg/cm 2 /s, (20) 



where t = 10 4 t4 s. This is consistent with the optical 
flux observed at early time, see figure Q] (note that our 
derivation holds only for r p h ~> R*, i.e. for t > 10 3 s). 

A more detailed study of the optical-UV early light 
could provide interesting constraints on the size and den- 
sity distribution near the surface of the progenitor star. 
Such a detailed analysis would require taking into ac- 
count effects neglected here (e.g., photon diffusion, which 
may be important on > 1 day time scale, anisotropy, etc.) 
and is beyond the scope of the current manuscript. 

Note that we have neglected the effects of photon diffu- 
sion in the above derivation. We do not expect diffusion 
to play an important role, due to the following argument. 
The size of a region around r(S m , t) over which diffusion 
has a significant effect is A(S m ,t) = y / ct/3Kp(5 m , t). 
Thus, the value of 6 m up to which diffusion affects 
the radiation field significantly, Sp, is determined by 
A(<5 m = S D ,t)/r(S m =5 D ,t) = 1, which gives 

r(6 m = 5 D ,t) 



= 1.2 f, 



-0.005 



p0.05 



\0.03 d 



-0.03 



(21) 



r ph . ' {M ej /M Q f 

This implies that photon diffusion is not expected to sig- 
nificantly modify the predicted light curve (Applying, 
e. g., the approxima te self-similar diffusion wave solutions 
of IChevalier] ri992. to our density and velocity profiles 
yields a luminosity that differs by < 25% from that de- 
rived neglecting diffusion). 

2.4. Radio emission 

For the massive wind discussed in our model, the syn- 
chrotron self-absorption optical depth is very large at 
radio frequencies. The characteristic frequency of syn- 
chrotron photons emitted by the lowest energy, j e ~ jt , 
electrons is v m ss j^eB /2Ttm e c, where the magnetic field 
is given by B 2 /8ir — espv 2 - Using eqs. (fl"0|) and (fTTj) we 
have 



v m « 5 x 10 u A- 2 e^ i a _ 1 c»_ 1 ^(m/10)- B / 6 t 



4/3 



-7/3 



Hz. 

(22) 

On a timescale of a few days we expect therefore the radio 
frequency, v ~ 10 GHz, to be in the range v m ^v<v c 
(see eq.©), i.e. we expect the Lorentz factor of elec- 
trons dominating the emission and absorption of radio 
waves to be in the range 7t < < 7c (where j c is 
the Lorentz factor of electrons with cooling time com- 
parable to the expansion time). In this case, the syn- 
chrotron self-absorption optical depth is given by r v = 
[e 3 B /2^ v (m e cv) 2 ]n e (^ v )A, where n e (7„) is the number 
density of electrons with Lorentz factor 7„ and A is the 
thickness of the shocked wind shell. For a power-law dis- 
tribution, dn e /c?7e °c 7~ 2 , we have n e (7„) « ti^t/Iv, 
and using n e A = pr/2m p for the electron column den- 
sity we have 



10 



~A 



e 1/3 (m/10) 5 / 3 ( - ) r 10/3 

1 W) V 10 GHz/ d ■ 
(23) 

This large optical depth leads to a strong suppression 
of the radio synchrotron flux, compared to the X-ray 
synchrotron flux given in eq. $7$ , 

{yL v )r&dio - „ a A 



(VL V \ 



i 10" 



x 



1/4 

e B.-l £ e -1 



E. 



-1/3 
fc,49 



(m/10) 



-11/12 



() 



The flux ratio measured at few to 20 days i s 
{vL v ) TaAio /{vL v )x » lO^ 2 5 (|Soderberg et alJ 12006ft . 
Thus, as noted in the introduction, the observed radio 
flux is higher than predicted by this model, and should 
be explained as a different component (see also §3.3[) . 

3. COMPARISON WITH OTHER AUTHORS 
3.1. Ghisellini et al. 2006 

iGhisellini etaTI (2006) argue that the shock breakout 
interpretation is not valid, since the optical-UV emission 
at few x 10 3 s is higher than the extrapolation to low fre- 
quencies of the Planck spectrum which fits the thermal 
X- ray emission at the sam e time. As we have pointed out 
in (jCampana et al. Il2006l ). and as explained in detail in 
§ [2j the thermal X-ray emission and the early optical- 
UV emission originate from different regions. The ther- 
mal X-rays originate from a (compressed) wind shell of 
mass ~ 10 _6 M Q , while the optical-UV emission origi- 
nates from the outer shells of the (expanding) star. At 
~ 10 4 s the optical-UV radiation is emitted from a shell 
at a "depth" of ~ 10 _4 M Q from the stellar edge. The 
optical-UV and the thermal X-rays need not correspond 
to the same Planck spe ctrum. As we have explained in 
(jCampana et al. 1l2006f) . and in more detail in §[2] if the 
explosion had been isotropic, the thermal X-ray emis- 
sion should have disappeared altogethe r on a time scale 
of few hundred seconds, at which time IGhisellini et al. I 
(2006) compare the X-ray and UV emission. The fact 
that the thermal X-ray emission is observed over a few 
thousand seconds can be accounted for by assuming an 
anisotrop ic explosion. This assu mption, made explicitly 
both in (ICampana et al. I I2006T ) and here, is ignored by 
IGhisellini et al. I (|2006f ). whose criticism is framed in the 
context of a spherical model, which we never assumed. 

Concerning their own model, in their preferred expla- 
nation for the properties of this burst IGhisellini et al. I 
(2006), unlike in our model, do not include any dynamics 
in their scenario, and adopt different ad-hoc parameters 
(e.g. for the electron distribution) at different times to 
account for the observations. 

3.2. Li 2006 

iLi I (|2006h argues that in order to explain the "tem- 
perature and the total energy of the blackbody compo- 
nent observed in GRB 060218 by the shock breakout, 
the progenitor WR star has to have an unrealistically 
large core radius ... larger than 100i? Q ". The results of 
ILi I (|2006[ ) are in fact consistent with our analysis. The 
"core radius", R c , adopted by this author refers not to 
the hydrostatic core radius of the star, but rather to the 
radius at which the optical depth equals 20. In non-LTE 
modeling of the winds of WR-stars, R c is located within 
the wind, typically near the sonic point. In fact, in some 
cases the wind velocity is already supersonic at R c (e.g. 
Hamann & Koesterke 1998). R c is much larger than the 
hydrostatic core radius Rhc obtained in evolutionary cal- 
culations (e.g. Schaerer & Maeder 1992). As shown by 
Nugis & Lamers (2002), R c /Rh c varies with the spectral 
type of WR stars, increasing from about 2 for early types 
to about 20 for late types (see table 5 of Nugis & Lamers 
2002). 

Since Li (2006) considers in his calculations models 
where the ratio between the wind photospheric radius 



R p h and R c is 1 < R p h/R c < 3 (using his notation, 
for b=5 the choice of e = 10~ 5 to 10~ 2 corresponds to 
1 < Rph/Rc < 3, see his fig. 3), his statement that 
R c ~ 100i?Q is required implies that a wind with a pho- 
tospheric radius of 7 x 10 12 cm to 2 x 10 13 cm is needed to 
account for the thermal X-ray emission of GRB 060218 
as shock breakout. This is consistent with our eq. Q. 

It is important to emphasize that the relevant radius 
is not R c , but rather the wind photospheric radius, R p h, 
where r = 1 . It is the photospheric radius which needs to 
be large in order to allow sufficiently large shock break- 
out energy. R p h is larger than R c by a factor of 2-10, 
depending on the model adopted for the wind velocity 
profile. Nugis & Lamers 2002 present late type models 
with R c /Rq > 30, and a large fraction of the stars ana- 
lyzed in Hamann & Koesterke (1998) have R c /Rq > 25 
(see their table 2). R p h larger than 100i? Q may therefore 
be obtained for stars with large mass loss rates. 

It should be emphasized that even if the wind photo- 
spheric radius required for GRB 060218 had turned out 
to be larger by a factor of a few than the largest R p h of 
known WR stars, this would not have been a strong ar- 
gument against the wind interpretation of GRB 060218, 
since GRB progenitors may well be more extreme than 
normal WR stars. Clearly, not all WR stars end their 
lives as GRBs. Moreover, it should be remembered that 
practically nothing is known about the mass loss on a 
day time scale prior to the explosion, which determines 
the wind density at the relevant radii. 

3.3. Fan et al. 2006 

iFan et ai~l (|2006l ) raised different criticisms in dif- 
ferent versions of their paper. In the first version, 
astro-ph/0604016vl, they argued that the presence of 
a massive wind would result in strong optical emis- 
sion at early time, which is inconsistent with observa- 
tions. While this claim was retracted in subsequent 
versions of their manuscript, it may be worthwhile to 
clarify this issue here. The synchrotron emission given 
in eq. (O holds for frequencies above the cooling fre- 
quency, v > v c . Since the cooling frequency is low 
for a massive wind, see equation l[5]). we expect in 
this model a similar X-ray and optical-UV luminosity, 
vf u «3x 10~ 13 (i/lday) _1 erg/cm 2 /s (corresponding to 
a flux of F uv « 1 x 10" 13 (t/lday)" 1 erg/cm 2 /s in the 
0.2/x band of Swift's UVOT). This predicted flux exceeds 
the observed flux at t < 10 4 s, hence the statement of 
I Fan et ai~l (|2006l ). However, as pointed out in § 12.21 
eq. ([7]) holds only for times t > t& ec ~ 10 4 s (and the lu- 
minosity is lower at earlier times). Moreover, at t < 10 s 
the large X-ray luminosity, Lx ~ 10 46 erg/s, would sup- 
press the synchrotron emission of the shock accelerated 
electrons. As explained in § 12.21 the cooling time of 
shock accelerated electrons due to synchrotron emission 
is short compared to the shock expansion time. At early 
times, t < 10 4 s, the ratio of the magnetic energy den- 
sity, esE/r 3 , to the X-ray energy density, Lx/4-7rr 2 c, is 
AiresEc/Lxr ~ 10~ 3 - 5 e Bi _ii? fc! 49(r/10 15 cm), which im- 
plies that the IC cooling time is much shorter than the 
synchrotron cooling time, suppressing the synchrotron 
emission in t he optical b and. 

Second. IFan et al. I (|2006h argue (see also 
iSoderberg et alJ 12006) that radio observations rule 
out the existence of a massive wind. As mentioned 
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in the introduction, indeed the radio and the X-ray 
afterglow are difficult to explain in the framework 
of a single shock. However, the radio observations 
alone cannot be used to rule out the existence of a 
massive wind, since radio observations do not allow 
one to determine the explosion parameters. This is 
illustrated by the fact that estimates for the kinetic 
energy based on modeling the radio data alone range 
from ~ 10 48 erg to ~ 10 50 erg, and that the ambient 
medium de nsity estimates range from ~ 10°cm~ 3 to 
~ 10 2 cm- ' 3 (iFan et al. 1 120061 ; ISoderberg et al.ll2006h . 

Finally, IFan et al. I p006) argue that in the presence 
of a massive wind, the radio flux would be higher than 
observed. As explained in § 12.41 our problem is quite the 
opposite: our model flux would be too low to account 
for the observed radio emission. ISoderberg et~aI1 (f200fih 
and IFan et ai~l (|2006h have chosen to construct mod- 
els of GRB060218 which concentrate on accounting for 
the radio emission, while attributing the X-ray afterglow 
to a continued activity of the source of an unexplained 
nature, and de-emphasizing the importance of prompt 
X-ray emission and the early optical-UV emission. How- 
ever, the radio emission repres ents a negligible fractio n 
of the total energy. Here, and in lCampana et al.~l (|2006|) . 
we have adopted a different approach, which is that the 
prompt (thermal and non-thermal) X-ray emission, the 
early optical-UV emission, and the late X-ray afterglow 
can all be explained within the context of the same 
model, which is based on the energetics of the early 
phases of the explosion. We argue that it is the radio 
afterglow, rather than all the other components, which 
plays a lesser role, and which may be attributed to a dif- 
ferent component. Given the very low energy of the radio 
emission and the the large anisotropy expected in the ex- 
plosion, such an additional low energy radio component 
is imaginable, whose role is unlikely to be important in 
determining the characteristics of the early high energy 
emission. 

4. DISCUSSION 

We have discussed a comprehensive model of 
the early X-ray and optical/UV behavior of the 
GRB060218/SN2006aj system, which provides the quan- 
titative justification fo r the interpretation outlined in 
iCampana et al. I ((2006), as well as a number of additional 
points. The most exciting features of this event were that 
it showed a strong thermal X-ray component as well as a 
strong optical/UV component in its early phases, at t < 1 
day, transitioning later to a more conventional X-ray and 
optical afterglow, and a radio afterglow. We have shown 
that the early X-ray/O/UV behavior can be understood 
in terms of a mildly relativistic radiation-mediated shock 
which breaks out of a (Thomson) optically thick wind 
produced by the progenitor star, leading to the observed 
thermal X-rays (see § 12. lj) . The early optical/UV be- 
havior arises as the shocked stellar envelope expands to 
larger radii (§ 12. 3[) , and the X-ray afterglow arises from 
synchrotron and inverse-Compton emission of electrons 
accelerated by the propagation of the shock further into 
the wind (§ 12. 2p . A detailed analysis of the optical-UV 
emission may therefore provide stringent constraints on 
the progenitor star. 

The thermal X-ray emission requires a mildly rela- 
tivistic shock, ~ 0.8, carrying > 10 49 erg, driven 



into a massive wind characterized by a mass loss rate 
of a few xl0~ 4 M Q /yr for a wind velocity of 10 3 km/s. 
The later X-ray afterglow is consistent with emission 
due to the propagation of this shock into the wind be- 
yond the breakout radius, where the shock becomes 
collisionless. This situation is very similar to that of 
GRB980425/SN1998bw, for which the X-ray and radio 
afterglow are interpreted as emission from a shock driven 
into a (much less massive) wind by the eje ction of a 
shell o f energy 10 49 7 erg and velocity — 0.8 (Waxman 
2004b], and references therein). 

The early optical-UV emission is consistent with the 
expansion of the outer part, SM ~ 10 _3 M Q , of the stel- 
lar envelope, which was heated to a high temperature 
by the radiation dominated shock as it accelerates in 
its propagation towards the stellar edge. It is impor- 
tant to note, however, that the acceleration of the shock 
near the stellar surface is not sufficient to account for 
the large energy, > 10 49 erg, deposited in the mildly 
relativistic, /3 ~ 0.8, component, which is required to 
have a mass of ~ 10 _5 Mq. For the parameters inferred 
from the SN2 006aj light curve, E ei ~ 2 x 10 51 erg and 
M e j ~ 2M (jMazzali et al.l [2006] ) . the energy predicted 
to by carried by > 0.8 is less than 10 45 erg (see fig. 6 of 
iTan. Matzner fc McKee1l200lD . This suggests that the 
mildly relativistic component is driven not (only) by the 
spherical SN shock propagating through the envelope, 
but possibly by a more relativistic component of the the 
explosion, e.g. a relativistic jet propagating through the 
star. 

An important factor in our interpretation of the early 
X-ray emission is the anisotropy of this shock, which 
leads to a characteristic timescale of the thermal X-ray 
emission controlled by the sideways pattern expansion 
speed, rather than by a simple radial line of sight ve- 
locity. A significant anisotropy of the mildly relativistic 
shock component is expected for various reasons, e.g. as 
a result of being driven by a jet, or due to propagation 
through an anisotropic envelope and into an anisotropic 
wind caused by fast progenitor rotation, etc. The sim- 
plest explanation for the prompt gamma-ray behavior 
may be that it is due to a relativistic jet, which can 
contribute to the anisotropy of the mildly relativistic 
shock propagating through the stellar envelope and the 
mildly relativistic ejected shell. In this case, the contri- 
bution of the relativistic jet to the late (t > 1 day) X- 
ray and optical aft erglow is not the dominant effect (c.f. 
iLiang et al. 1 12006). However, a possible explanation of 
the prompt non-thermal emission, both in GRB980425 
and GRB060218, may be that it is due to repeated 
IC sca tteri ng at breakout, a s suggested in IWang et al. I 
(|2006h (c.f. iDai et al. ir2006h . In this case there may be 
no highly relativistic jet- it may either be mildly rela- 
tivistic to begin with, or it may have been relativistic 
but it was choked, only its mildly relativistic bow shock 
emerging from the star. 

If our interpretation is correct, the importance of the 
early thermal X-ray component is that it represents the 
first detection ever of the breakout of a supernova shock 
from the effective photosphere of the progenitor star. 
Since this is a GRB-related supernova, a strong stel- 
lar wind can be expected, which results in the breakout 
photospheric radius being in the wind, rather than in the 



outer atmosphere of the star. This provides a potentially 
valuable tool for investigating the physical conditions, 
mass loss and composition of the long GRB progenitors 
in the last few days and hours of their evolution prior to 
their core collapse. 
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